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ABSTRACT 


Sound  fields  in  wedge-shaped  ocean  layers,  modeling  con- 
ditions on  the  continental  shelf,  have  been  studied  at  the 
Naval  Postgraduate  School  in  the  last  few  years  using  the 
method  of  images.  These  studies  are  carried  further  in  the 
present  work.  The  method  is  implemented  in  different  en- 
vironmental conditions.  This  thesis  examines  the  influence 
of  several  parameters  on  the  sound  field  for  downslope  pro- 
pagation in  a  wedge-shaped  fluid  of  speed  of  sound  C2  over- 
lying a  slow  bottom  of  speed  of  sound  c\ .  On  the  basis  of 
gualitative  and  semi-guant i tat i ve  analysis  of  the  behavior 
of  the  pressure-depth  profile  for  various  qeometrical  and 
physical  parameters,  we  can  conclude  that: 

1.  A  defined  distance,  the  "characteristic  distance" 
XQ  =  tt/(2)<2  s*n  9o  tan  &)'  where  cos  9Q  =  c^/^/ 
l<2  =  w/c2»  and  6  is  the  vertex  angle  of  the  wedae, 
has  physical  meaning  as  a  useful  scaling  distance. 

2.  The  distance  of  the  source  from  the  apex,  in  terms 
of  the  X0/  plays  a  major  role  in  determining  the 
downslope  sound  field. 
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I.   INTRODUCTION 

A.   SOUND  PROPAGATION  IN  SHALLOW  WATER  CHANNEL 

Experimental  investigations  of  sound  propagation  in 
shallow  water  channels  have  been  done  by  several  investi- 
gators. Shallow  water  propagation  is  of  interest  because  of 
the  applications  to  coastal  defense.  These  investigations 
are  expensive  and  time  consuming.  The  use  of  a  computer 
model  should  provide  a  relatively  inexpensive  alternative  to 
observation . 

One  of  the  technigues  uses  normal  mode  theory.  The 
normal  mode  theory,  which  was  introduced  and  developed  by  C. 
L.  Perkeris  [Ref.  1],  gave  the  exact  solution  in  water  of 
constant  depth.  Further  development  of  normal  mode  theory 
was  made  by  L.  Brekovskikh  [Ref.  2] ,  who  initiated  pressure 
as  an  integral  involving  Bessel  functions  and  solution  of 
the  normal  mode  eguation.  Another  theoretical  approach  to 
sound  propagation  in  a  horizontally  stratified  ocean  of  con- 
stant depth  is  given  by  the  method  of  multiple  scattering 
[Ref.  1],  With  this  method,  all  the  previous  theories  can 
be  simplified  by  conversion  into  an  asymptotic  form  which  is 
valid  when  the  acoustic  wavelength  is  small  compared  to  the 
distance  over  which  the  sound  speed  varies  appreciably. 
These  theories  agree  with  the  laboratory  experiments. 
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For  a  water  channel  with  a  small  bottom  slope,  the  sound 
field  may  be  expressed  approximately  in  terms  of  adiabatic 
normal  modes.  To  facilitate  prediction,  R.D.  Groves,  Anton 
Naql,  H.  Uberall,  and  G.  L.  Zauer  [Ref.  3]  modeled  a  wedge- 
shaped  isovelocity  ocean  with  a  linearly-sloping,  perfectly- 
rigid  ocean  floor  using  adiabatic  normal  modes.  For  a  pene- 
trable bottom  the  normal  node  description  fails  when  modes 
propagating  upslope  encounter  the  "critical  depth"  (Hc), 
defined  as  the  depth  where  the  associate  mode  changes  from 
fully  trapped  within  the  water  channel  to  radiating  energy 
into  the  bottom  (cut  off)  [Ref.  4-6] .  The  parabolic  egua- 
tion  can  be  used  to  explain  the  mechanism  of  sound  energy 
radiation  into  the  bottom  [Ref.  7,8] .  Such  an  equation  was 
studied  by  F.  B.  Jensen  and  W.  A.  Kuperman  [Ref.  9] ,  with 
predictions  that  satisfactorily  agreed  with  the  experimental 
results  for  small  ray  angles.  With  some  restrictions, 
normal  mode  theory  is  applicable  for  sound  propagation  in 
the  wedge-shaped  fluid  with  a  fast  bottom.  The  parabolic 
equation  is  good  for  fast  and  slow  bottom,  but  with  the 
restriction  that  horizontal  ray  angles  must  be  less  than 
20°. 

Another  technique  introduced  to  predict  the  propagation 
of  sound  in  the  wedge  is  the  method  of  images.  This  method 
was  derived  from  the  simplest  case;  a  monof requency  point 
source  in  a  homogeneous  ocean  with  parallel  boundaries.  The 
total  pressure  is  the  sum  of  an  infinite  number  of  spherical 
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waves  from  an  infinite  set  of  images.  The  restriction  of 
this  method  is  that  it  does  not  generalize  to  the  case  of 
inhomogeneous  media  or  non-olanar  boundaries.  In  this  work, 
this  method  will  be  studied. 

B.   THE  METHOD  OF  IMAGES 

In  1978,  Conpens,  Sanders,  Ioannou,  and  Kawamura  [Pef. 
10],  predicted  the  pressure  amplitude  and  phase  of  the  sound 
field  along  the  bottom  of  a  wedge-shaped  fluid  layer  of 
density  P]_  '  an^  speed  of  sound  c-j_,  overlying  a  fast  fluid 
bottom  of  density  Q2f  and  speed  of  sound  c2  >  c^  by  applying 
the  method  of  inaqes  in  a  computer  program  implementation. 
In  1934,  Baek  [Ref.  11],  and  LeSesne  [Ref.  12],  implemented 
further  improvements.  Baek  '  s  computer  program,  WEDGE,  and 
LeSesne 's  computer  prooram  XSLOPE  were  validated  for  several 
cases.  WEDGE  was  developed  for  two-dimensional  upslone  pro- 
oaoation  (the  source  and  received  are  in  the  same  vertical 
olane  pernendicular  to  the  shore  line,  and  the  receiver  is 
closer  to  the  apex  than  the  source  (Figure  1.1a))  and  down- 
slope  prorogation  (the  source  is  closer  to  the  apex  than  the 
receiver  (Figure  1.1b)).  XSLOPE  was  developed  for  upslope, 
downslope,  or  cross-slope  propaaation  (the  source,  receiver, 
and  apex,  are  not  necessarily  in  the  same  plane  perpendicu- 
lar to  the  shore  line  (Figure  1.1c)).  In  both  programs, 
Baek  and  LaSesne  assume  that  the  fluid  in  the  wedge  and 
fluid   in    the   bottom    have    constant    densities,   that 
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the  speed  of  sound   is  constant,  and  that  the   interface  be- 
tween the  fluids  and  the  surface  is  smooth. 

In  both  WEDGE  and  XSLOPE,  all  distances  are  scaled  in 
units  of  the  "dump  distance."  A  dump  distance  X,  as  stated 
in  Reference  10,  is  the  distance  from  the  apex  measured 
alona  the  interface  at  which  the  lowest  mode  attains  cutoff. 
If  the  wedge  angle  is  3  (Figure  1.1c),  then 


X  =  JL/2 (i  i) 

k]_  sin9c  tang  u,iJ 


'c 


=  arccos(c1/c2)  (1.2) 


where  k^  is  the  wave  number  in  the  wedge  and  9C  is  the  crit- 
ical grazing  angle  for  reflection  of  sound  from  the  bottom. 
For  3  <  <  1 

X  =  H  tan  6  (1.3) 

This  scaling  distance  negates  the  necessity  of  specifying 
frequency . 

C.   COMPUTER  PROGRAM  DSLOW 

At  the  start  of  the  work  reported  in  this  thesis,  a  com- 
puter program  was  obtained  [Ref .  13] ,  which  is  an  extension 
of  the  WEDGE  and  XSLOPE  for  downslope  configuration  with  a 
slow  bottom.  The  computer  model,  DSLOW,  developed  to  run  on 
a  desktop  computer  (Wang  2000),  uses  the  method  of  images  to 
predict    the    pressure    amplitude   and    phase    anywhere 
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within  the  wedge  fluid  overlying  a  slow  bottom  in  a  cross- 
slope  configuration.  A  geometrical  picture  of  this  configu- 
ration is  shown  in  Figure  1.1c. 

Mathematically,  the  model  used  in  WEDGE  and  XSLOPE  is 
applicable  in  any  condition.  But  consideration  must  be 
given  for  making  it  work  for  a  slow  bottom.  In  the  case  of 
a  fast  bottom,  the  dump  distance  has  a  physical  meaning. 
The  dump  distance  is  expressed  as  a  function  of  the  critical 
angle.  The  critical  angle  is  egual  to  arccos  (c^/c2).  In 
the  case  of  slow  bottom,  C]/c2  is  greater  than  1,  thus 
arccos  (c^/c2)  is  invalid;  therefore,  so  is  the  dump 
distance.  To  facilitate  the  scaling  factor,  a  "characteris- 
tic distance"  or  "scaling  distance"  is  introduced.  We  need 
the  scaling  distance  because,  with  this  distance,  our  model 
will  be  independent  of  freguency  as  in  the  fast  bottom  case. 
There  is  also  the  hope  that  the  use  of  a  scaling  distance 
will  allow  systematic  observation  of  the  pressure  field. 
This  scale  distance  XQ  is  the  distance  measured  along  the 
interface  from  the  apex  to  the  point  where  the  lowest  mode 
would  attain  cutoff  if  the  fluids  in  the  wedge  and  in  the 
bottom  were  to  be  interchanged.  The  characteristic  distance 
is  defined  by  the  following  eguation: 


tt/2 
k2  sin  8Q  tanB 


X0  =  " _.•_.   .__„  (1.4) 
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where  9o  =  arccos(c2/c1)  and  K2  =  w/c2  is  the  wave  number 
in  the  bottom. 

The  following  terms  will  be  used  throughout  (see  Figure 
1.1): 

8    =  wedoe  anale 

R^   =  distance  of  the  source  from  the   apex  in   units  of 


X 


o 


o 


=  distance  of  the  receiver  from  the  apex  in  units  of 
*o 

=  angle  of  elevation  of  the   source  above  the  bottom 

=  anole   of  elevation   of   the   receiver   above   the 
bottom 

=  distance  between  the  orojection  of  the   source  and 

receiver  on  the  shore  line,  scaled  by  X0 


p^/p2  is  the  ratio  between  the  density  of  the  fluid  in 
the  wedqe  ( p^ )  and  the  density  of  the  fluid  in  the 
bottom  ( p2 ) 

ci/c2  is  the   ratio  between   the  speed   of  sound  in  the 

wedqe  (c^)  and  the   speed  of   sound  in  the  bottom 

(cp)"    ^  fast  bottom  occurs  when   c2  >  Ci ;  a  slow 
bottom  occurs  when  c2  <  c^ 

The  purpose  of  this  research  is  the  following: 

1.  To  transfer,  test,  and  evaluate  DSLOW  program  on  the 
IBM  33  00; 

2.  To  obtain   numerical   and   Graphical   output   for   a 
number  of  cases;   and 

3.  To  attempt  to  develop  plausible  explanations  for  any 
significant  features  observed. 
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Figure  1.1   Geometry  of  the  wedge 
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II.  THEORY 

A.   GENERAL  VIEW  OF  A  WEDGE  PRESSURE  DISTRIBUTION  IN  THE 
DOWNSLOPE  CONFIGURATION 

A  general  picture  of  the  sound  energy  propagation  within 
the  wedge  in  downslope  direction  is  given  in  Figure  2.1.  If 
a  sound  source  is  placed  at  point  S,  ray  1  will  reach  the 
surface  at  point  P  with  an  incident  angle  a  with  respect  to 
the  normal  to  the  surface  at  this  point.  This  ray  is  re- 
flected by  the  surface  at  the  same  angle  but  with  the  phase 
180°  different.  (On  the  surface,  sound  pressure  is  zero 
everywhere.)  The  reflected  ray  reaches  the  bottom  with  an 
incident  angle  0  +  a  .  At  great  enough  distance,  ray  1 
never  reaches  the  bottom  again.  This  ray  does  not  contri- 
bute to  a  sound  pressure  field  at  the  bottom.  The  pressure 
at  the  bottom  should  be  very  small  according  to  the  ray 
theory  argument. 

Using  these  ray-tracing  methods,  an  estimated  profile  of 
the  pressure  amplitude  versus  the  receiver  depth  can  be 
made.  When  the  source  and  the  receiver  are  placed  near  the 
apex,  the  pressure  amplitude  is  zero  at  the  surface,  a  maxi- 
mum somewhere  within  the  wedge,  and  greater  than  zero  at  the 
bottom.  In  the  case  where  the  source  is  at  a  far  distance, 
the  pressure  amplitude  is  egual  to  zero  at  the  surface,  a 
maximum  somewhere  within  the  wedge,  and  zero  at  the  bottom. 
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Ray  tracing  will  only  give  a  rough  approximation,  not  an 
exact  solution,  but  ray  tracing  may  be  used  as  a  guide.  The 
method  of  images  calculates  the  exact  pressure  amplitude  at 
each  point  within  the  wedge  subject  only  to  the  assumption 
inherent  in  using  the  plane-wave  Rayleigh  reflection 
coefficients . 

B.   SOUND  PRESSURE  AT  A  POINT  IN  THE  WEDGE  DOWNSLOPE 
PREDICTED  BY  THE  METHOD  OF  IMAGES 

Let  the  source  be  a  scaled  distance  R\  from  the  apex  and 
at  an  angle  of  y  measured  from  the  bottom  of  the  wedge.  Let 
the  receiver  be  a  scaled  distance  R2  from  the  apex  and  at  an 
angle  6  measured  from  the  bottom. 

Using  Figure  2.2,  let  the  upper  half  family  of  images  be 
n  =  1,2,3,4,...  and  the  lower  half  family  be  n'  = 
1,2,3,4....  Calculating  the  field  resulting  from  source  and 
images  proceeds  along  the  lines  developed  in  [Ref.  14].  If 
(j>n  is  the  angle  formed  at  the  apex  between  the  nt^  image  of 
the  source  and  the  receiver,  then 

$1   =  2$  -  6  -  y 

$2  =  2  8  -  6  +  y 

4>3=  43-  6-  y 

<J>4=  46-  6+  y 


19 


Figure  2.1   Ray  tracing  in  the  wedge  downslope 
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or   more    generally 

<J>n   =     (n    +    l)S-6-Y  for    n    odd 

<j>n  =   nB  -  6  +  y      for  n  even 
Which  can  be  reduced  to: 

<j,n  =  {n+(l/2)  [l-(-l)n]  }B  +  (-l)ny  -  6         (2.1) 

or 

<frn  =  2  INT[   "  *  l       ]B  +  (-l)nY  -  5  (2.2) 

where  INT[  ]  denotes  the  largest  integer  which  is  egual  to, 
or  smaller  than  the  argument.  Using  the  same  method  for  the 
member  n'  of  the  lower  family  of  images  we  obtain: 

(j,n,  =  {n+(l/2)  [l-(-l)n]  }6  +  (-l)nY  +  5        (2.3) 
or 

4>n,  =  2  INT[   n  t    l       ]q    +    (-l)ny  +  6  (2.4) 

Using  the  geometry  of  Figures  2.2  and  2.4,  the  distance 
between  the  nt!l  and  n'*1"  images  to  the  receiver  is 
respectively 

rn  =  /  Rx2  +  R22  -  2R1R2cos$n  (2.5) 

and 

rn«  =  /R]/  +  R2^  -  2R1R2cos<}>n  (2.6) 
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Figure  2.2   Geometry  of  a  wedge  by  the  method  of  images 
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The  angles   9no   and   6n.0   for  the   ntn  and  n'th  imaqes 

respectively  are 

sin 


i_,0  =  arc  tan    - — 7- 

no  R2/P^ -  cos 


(2.7) 


and 


s  1  n  <j>n  1 


Li.  =  arc  tan  [  - — r=i 

no  R2/P]_  -  cos 


(2.8) 


Define  8nn?  and  9nin  as  the  angles  of  incidence  for  the 
mtn  bounces  from  the  bottom  for  the  n  and  n'  image 
respectively;  m  =  1,2,3....  (The  0th  bounce  is  the  last 
one  before  reaching  the  receiver.)  The  geometry  of  Figures 
2.2  and  2.3  oive  9nn  as  follows: 


21  =  920  " 

-  26  - 

5 

32  =  930  ' 

-  4Q    - 

5 

41  =  940  ' 

-    26  - 

6 

52  =  95U  " 

-  46  - 

6 

The  general  expression  is 
Qnm  =  9no  -  2m  6  -  6 
Using  the  same  method 


en*m 


•  n  1  0  -  2m  6  +  6 


The  maximum  number  of  bottom  bounces  of  the  ntn  and  nlth 


image    is 


fax 


=    M   =     INT     Un/26]    =     INT     Un,/26l 
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apex 


Figure  2.3   Geometric  development  of  reflection  angles 
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The    maximum    number    of    images    is 


rWv    =    N   =     INT     [it/3] 


'max 


The  reflection  coefficients   for  the  nth  and  n1  *-"  images 
for  a  plane  wave  are: 


R( 9nm ' 


Plcl 
P2C2 


_  ip 


nm 


Plcl 
P2C2 


+  i> 


nm 


(2.9 


and 


R(  9 


n'm} 


^2^2 


n  m 


Plcl  + 
P2C2 


^n' 


m 


(2.10) 


where 


and 


'nm 


/  1  -  (c1/c2)2cos2 


sin  9 


UUQ- 


nm 


(2.11) 


tn' 


n  m 


•  1  -  (c1/c2)2cos2  9 
sin  9n.m 


n '  m 


(2.12) 


The  contribution  from  the  upper  family  of  images  is 

N   .  M 

Pu  =  I  —  exp(-jkrn)(-l)INT[n+1)/2^  n  Rnm   (2.13 
n=l  rn  m=0 
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Figure  2.4   Geometry  of  symmetric  images 
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and  for  lower  family  of  images  is 


N   .  M 

Pu  =  X  — ,  exp(-jkrn.)(-l)INT[n+1>/2}  n  Rn-m 


n=l  rn 


(2.14 


m=0 


The  total  complex  pressure  is 


P(x)  =  Pu  +  PI 


(2.15) 
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III.   DSLOW  PROGRAM  IMPLEMENTATION 

A.   PROGRAMS  FEATURES 

Since  the  mainframe  graphics  computer  was  available,  the 
DISSPLA  graphical  program  was  used.  The  only  programming 
language  compatible  with  DISSPLA  is  FORTRAN.  The  numerical 
and  graphical  output  is  provided  by  this  program.  To  give 
the  pressure  amplitude  versus  received  angle  graphs,  two- 
dimensional  plotting  is  used. 

The  program  DSLOW  is  run  by  placing  the  point  source 
anywhere  in  the  wedge  and  then  placing  the  receiver  at  a 
distance  downslope  from  the  source.  The  receiver  position 
was  varied  from  zero  degrees  at  the  bottom  to  8  at  the 
surface.  High  resolution  plotting  was  achieved  by  dividing 
the  y-axis  (received  angle)  into  two  regions.  The  first 
region  covers  the  receiver  angles  from  zero  to  1/5  of  the 
wedge  angle.  In  this  region  A5  is  egual  to  8/100.  The 
second  region  covers  the  remaining  wedge  angle  with  A6  egual 
to  8/10.  This  method  provides  29  predictions  of  the  pres- 
sure amplitude.  Another  method  of  plotting  carried  out  was 
in  the  region  of  6  >  8/2,  A5  =  8/10,  and  in  the  region  of  5 
<  8/2,  A6  =  8/100.  This  method  provides  54  points  to  be 
plotted . 


28 


B.  NORMALIZATION 

The  main  goal  of  this  research  was  to  investigate  the 
profile  of  the  pressure  amplitude  as  a  function  of  a  number 
of  variables.  An  example  of  the  numerical  values  of  the 
pressure  amplitude,  the  normalized  pressure  amplitude,  and 
the  phase  at  each  receiver  position  is  displayed  in  Appendix 

C.  The  sound  pressure  becomes  smaller  as  the  receiver  is 
moved  away  from  the  source.  If  the  pressure  amplitude  were 
plotted  directly,  it  would  be  difficult  to  compare  the 
curves  at  near  distances  to  the  curves  at  far  distances 
since  at  the  near  distances  the  pressure  amplitude  is  much 
greater  than  the  pressure  amplitude  in  far  distance.  Thus, 
a  normalized  oressure  amplitude  is  needed.  The  normalized 
pressure  is  obtained  as  follows:   (see  Figure  3.1) 

We  know  that  the  sound  pressure  at  the  surface  is  zero 
and  that  the  sound  pressure  is  a  small  number  greater  than 
zero  at  a  point  near  the  surface.  The  first  non-zero  value 
of  pressure  P]_  is  at  the  receiver  angle,  6-j  =  93/10.  We  use 
this  first  calculated  non-zero  pressure  amplitude  as  the 
normalization  unit.   The  normalized  pressure  is 

PN  =  P( 5)/(P1)  (3.1) 

where  P(5)  is  the  pressure  at  any  point  within  the  wedge. 
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curve  of  the 
normalized  pressure  PN 
vs.  receiver  angle 


curve   of   the 

true   pressure    P    (S  ) 

vs.    receiver   angle 


pressure   amplitude 


Figure    3.1      Pressure    amplitude    normalization 
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C.  PROCEDURE 

Figures  3.2  through  3.7  represent  the  results  when  the 
receiver  distance  and  source  angle  are  fixed  and  the  source 
distance  and  receiver  angle  are  varied.  Figures  3.8  through 
3.13  reoresent  the  results  when  the  source  distance  and 
angle  are  fixed  and  the  received  distance  and  angle  are 
varied.  These  cases  will  be  the  foundation  of  our 
subseguent  discussions. 

The  solid  lines  indicate  the  fitted  curve  and  the  dots 
indicate  some  values  of  the  normalized  pressure  amplitude. 
In  DSLOW,  the  dot  appears  at  each  third  datum. 

D.  PROGRAM  IMPROVEMENT 

DSLOW  was  designed  to  provide  three-dimensional  graphs. 
For  example,  the  x-axis  represents  the  scaled  source 
distance,  the  y-axis  represents  the  scaled  received 
distance,  and  z-axis  represents  the  normalized  pressure 
amplitude.  To  simplify  the  presentation,  only  two- 
dimensional  graphs  were  presented  with  the  x-axis  the  nor- 
malized pressure  amplitude  and  the  y-axis  the  receiver  angle 
6.  All  curves  are  presented  with  the  data  fitted  with  a 
cubic  spline. 

The  DSLOW  program  was  executed  to  obtain  numerical 
results  of  the  phase  angle,  the  pressure  amplitude,  and  the 
normalized  pressure  amplitude  at  each  receiver  position. 
The  first  run  used  double  precision  for  accuracy. 
Difficulties   were   encountered  when   the  DISSPLA  subprogram 
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was  attached  for  making  the  graphical  output.  When  double 
precision  and  DISSPLA  were  not  successful,  the  single  pre- 
cision was  used,  resulting  in  round-off  error.  (See  Figure 
3.8  at  R2  =  10.0.  ) 
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Figure  3.2   Graphs  of  receiver  angle  5  versus  pressure 
amplitude  with  R2  fixed,  R^  varied 
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Figure  3.3   Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  Rj  varied 
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Figure  3.4   Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  R^  varied 
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Figure  3.5   Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  R^  varied 
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Figure  3.6   Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  R]^  varied 
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Figure  3.7   Graphs  of  receiver  angle  6  versus  pressure 


am 


plitude  with  R2  fixed,  Rj  varied 
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Figure  3.8   Graphs  of  receiver  angle  5  versus  pressure 
amplitude  with  R2  fixed,  R^  varied 
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Figure  3.9   Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  Rj  varied 
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Figure  3.10   Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  R^  varied 
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Figure  3.11   Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  F^  varied 
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Figure  3.12   Graphs  of  receiver  angle  6  versus  pressure 
amplitude  with  R2  fixed,  R±    varied 
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Figure  3.13 


Graphs  of  receiver  angle  6    versus  pressure 
amplitude  with  R2  fixed,  R^  varied 
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IV.   DISCUSSIONS 

\.       GRAPHICAL  OUTPUT 

The  graphs  of  normalized  pressure  amplitude  as  a  func- 
tion of  receiver  angle  5  were  investigated  for  various 
source  distances  R\  and  receiver  distances  R2  (Fig.  3.2  to 
3.13)  while  the  other  parameters  are  held  contstant.  For  a 
given  wedge  angle  3  and  sufficiently  small  source  distance 
(Fig.  3.8  and  3.12),  at  all  receiver  distances,  the  pressure 
increases  uniformly  towards  the  bottom.  For  greater  source 
distances,  (Fig.  3.9  -3.11)  the  pressure  attains  a  maximum 
within  the  wedge  for  all  receiver  distance. 

As  the  receiver  distance  is  increased  (Fig.  3.9),  a 
pressure  minimum  develops  between  the  maximum  and  the 
bottom.  An  important  property  of  the  curves  of  pressure 
versus  receiver  angle  when  there  is  a  maximum  and  minimum  is 
that,  at  a  specific  receiver  distance,  the  pressure  above 
the  minimum  can  be  extrapolated  to  zero  pressure  on  the 
bottom.  (See  Fig.  3.10  with  R2  =  32.)  This  receiver  dis- 
tance is  called  the  "transition  point."  So  far,  we  do  not 
know  the  properties  of  the  transition  point.  We  use  the 
transition  point  for  indicating  the  behavior  of  the  curves 
when  the  parameter  involved  is  varied.  The  transition  point 
appeared  twice  in  some  cases,  but  in  the  following 
discussions  the  first  transition  point  is  the  only   point  we 
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will  be  concerned  with.    (See  Fig.  4.1  for  transitions  cor- 
respond to  R2  =  4.6  and  6.4.) 

B.   GRAPHS  CLASSIFICATION 

The  development  of  curves  with  the  source  distance  (R^) 
and  the  receiver  distance  (R2)  as  variables  was  observed. 
As  R\  or  R2  are  varied  the  curve  changes  from  a  linear  curve 
to  a  curve  with  an  observable  minimum  (Fig.  3.9,  R2  =  5.0) 
and  finally  to  a  curve  without  a  minimum  (Fig.  3.9,  R2  = 
9.0).  Three  different  types  of  curves  resulted  from  the 
series  of  two-dimensional  plotting.  They  are  described 
below: 

1  .   Type  1  Curves 

Type  1  curves  (Fig.  4.2)  are  those  where  the  sound 
pressure  is  equal  to  zero  at  the  surface  and  maximum  at  the 
bottom  and  is  almost  linearly  dependent  on  depth.  These 
curves  are  most  pronounced  when  the  source  distance  is  much 
smaller  than  the  characteristic  distance.  The  closer  the 
source  is  to  the  characteristic  distance,  the  more  nonlinear 
the  curves  (Figs.  3.3  and  3.4). 

2 .   Type  2  Curves 

Type  2  curves  (Fig.  4.3)  are  those  where  the  sound 
pressure  is  zero  at  the  surface,  maximum  somewhere  between 
the  surface  and  the  bottom  with  no  minimum.  These  types  of 
curves  are  generated  when  the  source  is  placed  at  a  point 
much   greater   than   the  characteristic   distance.    Type   2 
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Figure  4.1   The  plots  where  there  are  two  transition  points, 
R2=4.6  is  the  first  and  R2=6.4  is  the  second 
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curves  indicate   that  the  sound   energy  in  the  wedge  is  well 
collimated  and  that  reflection  is  negligable. 
3.   Type  3  Curves 

Type  3  curves  (Fig.  4.4)  are  those  that  have  a  min- 
imum pressure.  These  curves  occur  when  the  source  is  a  dis- 
tance slightly  greater  than,  or  less  than,  the  charac- 
teristic distance.  Tables  1,  2,  and  3  of  Appendix  D  show 
the  receiver  positions  at  the  first  transition  points. 
Three  different  values  of  8,  two  different  values  of  pi/P2» 
and  two  different  values  of  ^\/^2  were  used  in  making  these 
tables.  The  transition  point  did  not  occur  when  8  =  15°, 
p-L/P2  =  0.90,  c^/c2  =  1.10.  An  explanation  can  be  offered 
using  the  fact  that  for  these  particular  sound-speed  and 
density  ratios  an  angle  of  intromission  exists  [Ref.  15]. 
Since  the  angle  of  intromission  is  the  grazing  angle  at 
which  the  sound  energy  is  completely  transmitted  into  the 
slow  bottom,  it  is  plausible  that  no  transition  point 
occurs . 

C.   TRANSITION  POINT 

By  varying  the  wedge  angle  8  in  small  increments  A8  = 
0.5°  starting  with  8  =  5°,  and  ending  at  8  =  7°,  it  was 
found  that  transition  occurs  for  source  distances  within  the 
range  from  1.0  to  1.5. 

Ror  R;l  <  1.0,  no  transition  point  was  observed;  the 
curves  are  the  Type  1.   For  1.0  <  Ri  <  1.5,  the  evolution  of 
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curves  as  the  receiver  distance  varied  can  be  explained  as 
follows:  first,  the  receiver  is  placed  near  the  source  and 
gradually  it  is  shifted  further  from  the  source.  The  mini- 
mum in  the  pressure  decreases  reaching  the  point  where  the 
curves  extrapolate  to  zero  (the  first  transition  point). 
Further  detailed  observations  were  made  on  this  particular 
facet  by  varying  the  source  distance  and  the  receiver 
distances.  The  results  of  these  observat i<)  is  are  tabulated 
ere  < raphed  in  Appendices  D  and  E.  When  the  receiver  is 
moved  away  from  the  source,  the  minimum  will  reach  a  minimum 
pressure  then  the  pressure  increases  until  it  reaches  the 
point  where  the  curves  again  can  be  extrapolated  to  zero, 
this  is  the  second  transition  point  (See  Fig.  4.1). 

For  R\  >  1.5,  there  will  be  no  transition  point.  The 
curves  are  the  Type  2. 

The  transition  point  as  a  function  of  source  angle  can 
be  observed  using  the  tables  in  Appendix  D.  In  most  cases 
the  greater  the  source  angle  -y,  the  closer  the  transition 
point  is  to  the  apex.  Graphs  of  transition  point  as  a  func- 
tion of  R^  (Appendix  E)  indicate  that  the  smaller  3  the  more 
regular  the  curves.  This  is  easy  to  understand  because  the 
smaller  8,  the  more  accurate  the  observation  of  transition 
point;   the    greater   3    the   less    accurate    the   data. 
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Figure  4.2   Type  2  curves,  indicating  a  pressure 
amplitude  nearly  linear  with  depth 
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Figure  4.3   Type  2  curves,  indicating  a  well-collimated 
sound  field  as  the  source  away  from  the  apex 
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Figure  4.4   Type  3  curves,  indicates  the  presence  of 
reflection  and  refraction  near  the  bottom 
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D.   PARAMETER  VARIATIONS 

Variation  of  parameters  was  done  by  changing  one 
parameter  of  interest  while  all  others  were  held  constant, 
for  fixed  source  and  receiver  distances,  and  plotting  the 
receiver  angle  versus  normalized  pressure  amplitude. 

The  parameters  6,  P]/P2'  and  C\/Q2'  were  held  constant 
and  the  pressure  amplitude  was  plotted  for  various  6,  R^  , 
and  R2 •  The  y's  are  set  at  6/4,  6/2,  and  3  6/4.  Variations 
in  the  shore  distance  (YQ)  can  be  made  because  the  program 
is  available,  but  to  simplify  the  investigation,  Y0  was  set 
egual  to  zero  for  all  plots  (Fig  4.2  is  included  as  an 
example  for  Y0  *    0). 

1.  Variations  of  6 

Initially,  the  values  of  6  investigated  were:  6°, 
10°,  and  15°.  The  major  effect  created  by  altering  the 
value  of  6  is  that,  for  the  same  values  of  R\ ,  p^/s>2'  and 
c^/c2f  the  smaller  6,  the  shorter  the  transition  point  (see 
Figs.  3.9  and  3.10,  Tables  1,  2,  and  3  of  Appendix  D) . 

2.  Variations  of  y 

The  variations  of  y  from  y  =  6/4  (the  source  is 
placed  near  the  bottom)  to  y  =  36/4  (the  source  is  placed 
near  the  surface)  are  presented  in  Tables  1,  2,  and  3  indi- 
cated that  the  greater  y  the  shorter  the  transition  point. 
It  is  not  always  true,  for  instance  in  Table  1  at  6  =  6°,  R]_ 
=  1.50,  the  greater   y  the  longer  transition  point,   for  the 
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rises.   (See  Appendix  D. ) 

3.   Variations  of  c^/c2  and  p^/p2 

Variations  of  the  acoustical  parameters  ci/°2  an<^ 
pi/p?  w^re  done,  but  did  not  give  a  significant  variation  of 
the  sound  pressure  profile. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

The  program  DSLOW  gives  good  plots  representative  of  the 
sound  energy  distribution  within  the  wedge.  The  sound  energy 
can  be  well  collimated  by  the  wedge.  This  phenomenon  is 
strongly  affected  by  the  source  position.  At  a  source 
position  close  enough  to  the  aoex,  sound  energy  is  distri- 
buted linearly  with  respect  to  the  depth.  As  the  source 
moved  away  from  the  apex,  the  distribution  of  sound  energy 
becomes  more  complex.  Sometimes  a  minimum  is  found;  this 
minimum  may  be  caused  by  the  presence  of  sound  energy 
reflected  by  the  bottom. 

The  source  position  plays  a  major  role  in  forming  the 
pressure  distribution  profile.  The  pressure  distribution  is 
also  very  sensitive  to  the  parameter  variation  at  small 
source  distances,  but  it  becomes  insensitive  at  large  source 
distance.  The  characteristic  distance  must  have  physical 
meanings  rather  than  just  an  arbitrary  number,  because  when 
the  source  distance  in  proximity  to  the  characteristic 
distance,  the  model  is  most  sensitive. 

The  model  is  restricted  when  the  single  precision  mode 
generates  round-off  error  and  rough  curves  which  do  not 
allow  for  accurate  analysis. 
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B.   RECOMMENDATIONS 

1.  Sincjle  precision  produces  aood  results,  but  failed 
in  some  cases.  Double  precision  would  improve  the 
program,  but  increase  the  execution  tine.  This  must 
be  done  by  running  the  program  in  double  precision, 
and  accumulating  the  result  in  sincjle  precision 
before  plottinq  the  data  by  DISSPLA. 

2.  It  is  sugoested  that  the  program  be  run  using  more 
realistic  parameters  and  observing  the  effects  on 
the  characteristic  distance  and  transition  point. 

3.  Further  study  validating  DSLOW  in  comparison  with 
exoerimental  results  is  suggested. 
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APPENDIX  A 
DSLOW  ALGORITHM 

The  pressure  amplitude  calculation 
N-!  =  INT[180,'P]  (eqn  A.l) 

KiX  = (eqnA.2) 

2tanp  tan[arccos(c-).Cj)] 

AL  =  a  K2  =  0.0001      (constant)  (eqn  A. 3) 

D2  =  Y02+  Rj2  +  R22  (eqn  A.4) 

R3  =  2RjR2  (eqnA.5) 


R8=V°2'  R3  cos[(N-l)P  +  Y  -6]  (eqn  A. 6) 


R9=VD2- R3  cos[(N-l)p  +  Y+  6]  (eqn  A. 7) 

S2  =  (-l)INT(Nl'2)  (eqn  A. 8) 

W1  =  2AL(c1,c2)2  (eqn  A. 9) 

SI  =  |{Rj  sin[(N-l)p  +  y]-2[INT((N-l)/2)p  (eqn  A. 10) 
+  R2sin[2INT((N-l),2)P-5]}|/R8 

CI  =  V/(1-SI2)  (eqn  A.l  1) 

1=81/0!  (eqn  A.  12) 

Wo  =  (-c2  +  cl'c2>  (eqnA.13) 
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Y=VW02  +  Wj^  (eqnA.14) 

Z=|W0|  (eqnA.15) 


Yj-^Y  +  Wq),^  (eqnA.16) 


Y2  =  V/(Y-W0),2  (eqnA.17) 

T  -Y 

Z,= ^ r  (eqnA.lS) 

1       (T-Y2r  +  Y^ 


Y, 

Z^  = ri r  (eqnA.19) 

-      (T-Y^  +  Y^ 

(T2.Y,)2-Y2 

(TM2r+Yi 

2Y,T 
Z^=  (T^-Y^V  (CqnA-21) 

N 

1 

Pi-  I  (-1)INT(V2)(Z5  cos(R8nKjX)  +  Z6  sin(R8nK1X},  R8n  (eqn  A.22) 

n=l 

N 
1 

P2=  I  (-DINT(N;2){-Z5  sin(R8nK1X)  +  Z6  cosfRg^XJ/Rgn  (eqn  A.23) 

n=l 

N 
1 

p3=  I  (-DINT{N;2){Z5  cos(R9nK1X)  +  Z6  sinfR^KjXJ/Rgn  (eqn  A.24) 

n=l 
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N 

1 

P4=  I  (-1)INT(X;2>{-Z5  sin(R9nK1X)  +  Z6  cos(R9nKjX}  R9n  (eqn  A.25) 

n=l 

P5=  Pl  +  P2  (eqn  A.26) 

P6=  P3  +  P4  (eqn  A.27) 


Ptot=RlVP5^P6Z  (eqnA.28) 
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APPENDIX  B 
DSLOW  PROGRAM 

******** ********************************************* ******************* 

**  THIS  PROGRAM  CALLED  DSLOW  IS  CALCULATING  THE  SOUND  PRESSURE  ** 
**  WITHIN  THE  WEDGE  OVERLYING  SLOW  BOTTOM  FLUID  AND  DOWN-SLOPE  ** 
**  DIRECTION  ** 

*********************************x*x ************************************ 

-     INTEGER    A,I,I1.M.N,S1,S2.N1.J,K,P 

REAL*4  6, Cc!C2j D,6l,62  GJPI .PI .P2-Ql-RliR2.T- 

******************** ***************************************************** 

**     INPUT  PARAMETERS  ** 

************************************************************************* 

C  B  =  WEDGE  ANGLE  (DEG) 

C  G  =  SOURCE  ANGLE  TDEG) 

C  D  =  RECEIVER  ANGLE  TDEG) 

C  Nl=  #  OF  IMAGE  POINTS 

C  Rl  =  SOURCE  DISTANCE  (IN  CHARACTERISTIC  DISTANCES)  , 

C  R2  =  RECEIVER  DISTANCE- (IN  CHARACTERISTIC  DISTANCES) 

C  YO  =  APEX  DISTANCE  (IN  CHARACTERISTIC  DISTANCES) 

C  Dl  =  RHO  1/RHO  2 

C  CC  =  C-l/C-2 

C  AL  =  ALPHA/K2 

C  A  =  #  OF  RECEIVER  POSITIONS 

********************************  **************************************** 

**  TMTTTAT      TMPIIT    K,    BFT  1TTDM";  ** 


INITIAL  JINPUT.&.REUTION 

r************************ 


******  **  *  ***  *  *  **  *************  ************************************  ******* 


C 
C 

c 

c 
c 

C 

c 

c 
c 

-      CC  =  l.ll 

AL_=  0.0001 

****** ******************************************************************* 

********************J^*K************************************************ 

K   =  0  , 

Nl  =  INT(180./B) 

T6  =  180. /PI 

B  =  B/T6 

G  =  G/tS 

C2  =  1:0**2 

=  AYg*Yp)+(Rl*Rl)+(R2*R2) 

~!*TAN(ACOS(l/CC)  )*TAN(B) ) 


YO 

=  0 

.00 

R2 

=  9 

.00 

DO 

111 

P=l,6 

B  = 

'  10 

.0 

B  = 

=  10 

.0 

G  = 

:  3;b/4 

G  = 

:  2*B/4 

G  = 

>  l*B/4 

D  = 

■■   0.0 

Rl 

=  0. 

.95 

Rl 

=  1 

.10 

Rl 

=  1 

:g 

Rl 

=  0 

R2 

=  1 

Ah 

R2 

=  1 

R2 

=  1 

.50 

R2 

=  2. 

.00 

R2 

=  D 

.00 

YO 

=  0 

.00 

Dl 

=  0 

.80 

AB, 

11  =  ACOSfl/CC) 

"2  =  TAN 7 TOO  IT 

3  =  2.*TQq2*'tQQ 

=  PI/TQQ3 
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D2  =  Y0*Y0+R?*R1?R2^R2 

01  =  1/DSQRT(2.0D00) 
C240   FORMATS  NOP    ' , 7X, ' THETA(N)  ' ,7X,'IMGE  SR  R8   ',7X, 
C    *     ' IMGE  SR  R9  ' ) 


C     WRITE  (6,240 
-241   FORMAT?' 

WRITE  (6,241 
-250   FORMAT  '  WEDGE  ANGLE  =    ',F5.2,1X,'  SOURCE  ANGLE  =   \F5.2) 

WRITE  (6 ,250]  B*T6,G*T6 
-251  TORMAT('  SOURCE  DISTANCE= ', F4 . 2, IX, ■ RECEIVER  DISTANCE=  ',F6.2,1X, 

*  ,'  SHORE  DISTANCE=  ',F4.2) 
WRITE  (6,251)  R1.R2.Y0 

-270   FORMAT(r  RHOi/RH62='  F5.2,5X, '  Cl/C2= ■ , F5 .2 , 5X, '  ALPHA/K2= ■ , F8.4) 

WRITE   6,270)  D1,CC,AL    '''''-' 
-271   FORMATS  K1XV   =   \F5.2) 

WRITE  (6,271)  T4 
-272   FORMATS '   ')' 

WRITE  (6.272) 
-800   FORMAT ( IX,' REC.POS'2X'  REC. ANGLE  ' ,2X, 'PRES. AMPLITUDE  '  ,3X, 

*  'PHASE  ANGLE  ' ,4X, 'NORM. PRESS ' ) 
WRITE  (6,800) 

-801  FORMAT {'  ') 
WRITE  6,801) 
DO  10   M  =  0,A 

D  =  M*B/A 

DD=B/A  , 

V  =2*B/10 


IF(D.Lt.V)  GOTO  110 


IF(D.GE.V)  GOTO  120 
•110     DX  =  D 

DO  15   J  =  1,1CL  ,, 

D  =  DXt(J;l)*B/(10*A) 

DD  =  B/ilO*A)      . 
120     SI  =  1.0 

DO  20   N  =  1,N1 


IF(Sl.GT.O)  T1(N)=(N-1)*B+G 
IF(Sl.LT.O)  Tl  N)=N*B-G 
SI  =  -  SI 

R8(N)  =  SQRT(D2-R3*C0S(T1(N 

R9{N)  =  SQRT(D2-R3*C0S(T1(N 

FORMAT(3XlI2  *5X;F5.2.4X.F6 .4 ,6 

WRITE  (6,310)  N,Ti(Nj,R8(N),R9(N 


C310 

C 

-20   CONTINUE 

PI  =  0.0 
P2  =  0.0 


DO  30  N  =  1,N1   ,  ,  vx 

Wl  =  2*C2*AL  . 
II  =  INT((N-l)/2) 
DO  40v  I  =  1.11 
S(I   = 


)  =  ABS[R1*SIN(T1(N)-2*I*B) 
+R2*SINi^*I*B-DJ  /R8(N) 
CI   =TSQRTU-(S(I)^ST1 ))) 
T  =  S(I1/D1 

wo  =  l-c2+tcm*c(im  V1 
y  =  s6rt7(wo^w6)+{w1*wi)) 

Z  =ABSiW0) 
IF(Y.LE,ZJ  Y=  Z  , 
Yl  =  Ql*s6RT(Y+W0h 
Y2  =  -Ql*SQRT(Y-w6) 
Zl  =  T-Y2 
Z2  =  -Yl   , 
:i2=  CMPLX(Z1R,Z2I)  . 


Z3  =  Zl/7zi*Zl4Z2^. 
Z4  =-Z2/(Zl*Zl+Z2*Z; 
Zl  =  T+Y2 
Z2  =  Yl 


Z13=  CMPLX(Z3R.Z4I) 

Z5  =  Z1*Z3-Z2*Z4 

Z6  =  Z1*Z4+Z2*Z3   v 
C  214=  CMPLX(Z5R,Z6I) 

E(I      =   Z5 

F(I)    =  Z6 
C400   ^FORMATf' ;nOP      '    5X,'I    ' , 5X, ' ARSIN(S (I )*T6) ' , 5X, ' E(I) ' , 

C         *WRITE    (^DSgftT(E?t)*E(I)+F(I)*F(I))'/5X,'T6*ATAN(F(I)/E(I))') 
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C420      FORMAT(3X.I2<6X,I2,5X.F6.4.10X,F6.4.3X,F6.4,3X.F6.4,10X,F7.4) 

C  WRITE    (6,420)    Nil.    ASlNTsTtWte.EliLHlL 

-      40    CONTINUE  S$Rt(ETl)*EUMl  J*Mf  J)  ;ATAN(F(I)/E(I)  ) 

21-0 

Z2   =   0 
Z3   =   0 

24  =   0 

25  =   1 
Z6   =   0 

IF(N,LE.2.00)    GOTO   50 
DO  45      I   =   1,    II 
Zl   =  E(I) 


U 


Z2  =  Pi 
Z3  =  Zl 
24  =  Z6X 

Z5   =  Z1*Z3-Z2*Z4 
Z6   =  Z1*Z4+Z2*Z3 
CONTINUE 
Zl 
Z2 
T 

Z3 
24 

Z6  =  Zl*Z4+Z2*Z:y 
PI  =  P1+S2*Z5/R8(N 
P2  =  P2+S2*Z6/R8(" 


m 

,  'RI 


C500  FORMAT ('    NO, OF    I.P    ' , 3X    ' RE(REFL)=    ' ,3X, ' IM(REFL)=    ') 

C  WRITE    (6y506) 

C510  F0RMAT(3X,I2<6X,F6,4.6X.F6.4) 

C  WRITE    (6,510}   N  S2*Zl.S2*Z2    ' 

C600  FORMAT (  '     I'  ,5X.  'ARCSIN(S(I))*T6='  ,3X,  'E(I)'  ,5X,  'F(I)'  ,5X,  'EF'  , 

C  *  5X    '    AN2    ' T 

C 


WRITE    (6.660) 

ii=n+i 


:(i)  =  :  drt'C 


S(I)'=   ABS(R1*SIN(T1(N)-2*(I-1)*B}      „ 
+,R2*S|N  C|*«nJ*B+0})/R9(N) 


T  ■  S(IT7DI 

WO   =   -C2+ctl)*C(I) 

Y   =   SQRT(TW0*W0)+(W1*W1)) 

Z   =   ABS(WQ) 

IF (Y.LE.Z)Y   =   Z    , 

Yl    =   Q1*SQRT(Y+W0h 

Y2  =  -Q1KSQRT(Y-W6) 

Zl  =  T-Y2  * 

Z2  =  -Y 


Z3  =  Zl7(Zl*Zl+Z2*Z2) 
Z4  =  -  Z2/(Z1*Z1+Z2*Z2 
Zl  =  T+Y2 


=  Yl 

=  Z1*Z3-Z2*Z4 

=  Z1*Z4+Z2*Z3 


,1)  =  Z6,  .  t%i 

AN1=  ASIN(S(I))*T6X   ,  _  ,  %. 

EF  =  SORT  E(I  *E(I)+F(I)*F  I)) 

AN2=  ATAN(F  I  /E (I)  *T6 
C610   FORMAT(2X.I2<5X/F6.4.5X<F6,4<5X.F6.4,5X,F6.4,5X,F7.4) 
C     WRITE  (6,610)  I,ANl,E(l),F(l  ,EF,AN2 
-   60  CONTINUE 


DO  80 


80 
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C700 

c 

C710 

C 

-   30 


C 
-810 


2  =  26^ 

=  T4*R9(N) 

3  =  COS?Tl 

4  =  -SIN(T) 

5  =  Z1*Z3-Z2*Z4 

6  =  Zl*Z4tZ2*Z3,  v 

1  =  P1+S2*Z5/R9(N) 

2  =  P2+S2*Z6/R9(N} 
LOWER  PATH  N0=>,3X, 'RE(REFL)=  ' , 3X , ' IM(REFL)=  ') 

x:i2rl2X,F6.4,12X,F6.4) 
,710)  N,$2*Z1,S2*22 


•10 


-811 

-31 

C 

-111 

***** 

** 

***** 

C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


FORMAT 

WRITE 

FORMAT 

WRITE 

CONTIN 

K  =  K+l 

DZ(K)=D*T6 

PZCK)=SQRT(Pl*PltP2*P2)*Rl 

AT(K  )=ATAN(P2/P1 ) 

WRITE  (6.810)  K,DZ(K),PZ(K) . ATAN(P2/P1 ) ,PN(K) 

FORMAT (6X, 13  11X,F5.2  12X, F7 .4 , 12X,F7 .4  10X,F7.4 

CONTINUE 

CONTINUE 

DO  31,L=1,29 

XP=PZ  23)  x  , 


)  31 ,L=1 

>=PZ  23)  , 

J(L)=PZ(L 
WRITE  (6,8liy"L,DZ(L).PZ(L),AT(L) .PN(L) 
FORMAT(2X,I3  7X  F5.2  8X,F9.6,8X,F9.6  5X,F9.6) 


iV) 


CONTINUE 

Y0=Y0+0.5 

R2=R2*0.50 

CONTINUE 

STOP 

END 

***************************************************************** 
***************************************************************** 

CALL  MEDBUF 

CALL  TEK618 

CALL  COMPRS 

CALL  NOBRDR 

CALL  PAGE[1! 

CALL  AREA2DI 

CALL  HEIGHTS . 

CALL  XNAME V NORMALIZED  PRESSURE  AMPLITUDES  \  29) 

CALL  YNAMEC' RECEIVER  ANGLE (DEG)$ ', 19) 

CALL  YTICKS(5V 

CALL  XTICKS(5 

CALL  GRAF  CO 

CALL  GRAF(0 

CALL  GRAF(0 

CALL  DOT 

CALL  GRID(2 

CALL  HEADINi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGI 

CALL  MESSAGi 

CALL  MESSAGI 

CALL  MESSAGI 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGI 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGI 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 

CALL  MESSAGi 


100,1.8,1) 


1  ,100 

;ioo 


,0.10.0.50.0. .5.0,15.0) 

,10. ,$o.o,0. ,5.0,10.6}' 

, 3. 0;i5.0;0.;3. 0,15.0) 

'REC. ANGLE  VS.  PRESSURES ', - 
'WEDGE  ANGLE:  $'100.8. .71) 

1 15. os »  ioo, 'abut ' , 'Abut ' ' 

'RH01/RH02:  $',100,8. ,6. 5 

'0.80  s'.ioo. 'Abut' , 'Abut 

'Cl7C2:  S'  160.8. , 6.0T 
'1.505' ,106. 'ABUT', 'ABUT') 
'SOURCE  ANGLE:  $' .100, 8. .5.5) 
'ABUT' , 'ABUT1} 
'ABUT' , 'ABUT' ) 
'ABUT' , 'ABUT' J 
'ABUT' ' ABUT' )v 
:  S'  I6O.8.  ,5  ) 

'ABuf ' , 'Abut ' 

, 'ABUT' , 'ABUT' 
. 'ABUT' , 'ABUT' ) 
,106, 'ABUT', 'ABUT* , 
.DIST .INtT:  $' .100.8. .4.5) 
' ,100, 'ABUT' , 'ABUT'? 
1 ,100, 'ABUT' , 'ABUT' 
1 ,100, 'ABUT' , 'ABUT' 
■100. 'ABUT'  'ABUT' , 
St.INCRMT  :  $',100,$. ,4.0) 

'ioo, 'abut' , 'Abut  Y 

,100^ 'ABUT' , 'ABUT' ) 
$T.INIT:  S' .100.8.  3.5) 
T,  100, 'ABUT' , 'ABUT') 


'11.25 

1  3.75$',100 
1  7.50S':i00 
1  SOURCE. DIST 
'0.25  S' ,100 
'1.10  $',100 
1.6Q$Tfi00 


'0.9 

'REC 

■0.5 

'1.0 

■3.00 

'3.00 

'REC.D 

'TWICE 

•2.00$ 

■SHO.D 

'1.00 
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C  CALL  MESSAGE 'SHO.DIST:  $ ' . 100 .8. ,3 .5) 

C  CALL  PARA3 

C  CALL  POLY3 

C  CALL  NOCHEK 

C  CALL  CURVE (PN,DZ, 29,0) 

Clll  CONTINUE   V 

C  CALL  ENDPL(O) 

C  CALL  DONEPL 

C  STOP 

C  END 


**  For  plotting  the  calculations  result  on  TEK618  or  Sherpa,  erase   ** 
**  the  proper  C's  in  DISSPLA  subprogram  above.  Put  C's  in  front  of   ** 
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APPENDIX  C 
NUMERICAL  RESULTS  OF  DSLOVV 


WEDGE  ANGLE  =    10. 0( 
SOURCE  DISTANCED. 75 
RK01/RH02=  0.80 
K1X   =   19.44 


SOURCE  ANGLE  =    2.50 
RECEIVER  DISTANCE=    9.00 
Cl/C2=  1.10 


SHORE  DISTANCE=  0.00 
ALPHA/K2=   0.0001 


REC.POS    REC. ANGLE    PRES .AMPLITUDE 


.0 
_.1_ 
0.20 
0.30 


0 

0 

0 

0 

80 

90 

00 

10 

20 

30 

40 

50 

60 

70 

18 

00 

00 

00 

00 

00 

7.00 

8.00 

9.00 

10.00 


0 
0 
0 
0 
0 
0. 

1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 

2. 

3. 


M 


.0071 
. .00739( 
0.007175 
0.006965 
0.006763 
0.006563 
0.006374 
0.006191 
0.006014 
0.005844 
0.00568 
0.00552 
0.005379 
0.005239 
0.005107 
0.004983 
0.004866 
0.004759 
0.004658 
0.004566 
0.004482 
0.004003 
0.003919 
0.003841 
0.003550 
0.002981 
0.002149 
0.001121 
0.000000 


PHASE  ANGLE 

NORM. PRESS 

-1.129606 

6.792701 
6.595129 

-1.108071 

-1.085083 

6.403362 

-1.061422 

6.215911 

-1.037255 

6.035660 

-1.011464 

5.857465 

-0.984591 

5.688235 

-0.957120 
-0.928227 

5.524890 

5.367038 
5.215738 

-0.898359 

-0.867409 

5.070300 

-0.835250 

4.931440 

-0.802517 

4.800441 

-0.768184 

4.675240 

-0.732669 

4.557538 

-0.696848 

4.446329 

-0.659586 

4.342858 

-0.622856 

4.247226 

-0.583461 

4.156956 

-0.545011 

4.074470 

-0.505302 

3.999775 

-0.113330 

3.572277 

0.201979 

3.497443 
3.427753 

0.413213 

0.547851 
0.632015 

3.168519 

2.660608 

0.680885 

1.917873 

0.697709 

1.000000 

-0.048596 

0.000029 
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WEDGE  ANGLE  =    10.00 
SOURCE  DISTANCED. ^5 
RH01/RH02=  0.80 
K1X   =   19.44 


REC.POS 


6 
7 
8 

J 
i 

13 

if 

17 
18 
19 
20 
21 
22 
23 

I! 

26 
27 
23 
29 


RE C. ANGLE 
.00 


SOURCE  ANGLE  =    2.50 
RECEIVER  DISTANCE=   4.50 
Cl/C2=  1.10 


SHORE  DISTANCE*  0.00 
ALPHA/K2=   0.0001 


PRES. AMPLITUDE 


.03?j 


.32 

.037361 
0.036906 
0.036460 
0.036019 
0.035594 
0.035178 
0.034766 
0.034371 
0.033981 
Q  .033602 
0.033229 
0.032367 
0.032511 
0.032163 
".031823 
.031492 
.031166 
.030346 
.030532 
.030226 
.027368 
0.024629 
0.021663 
0.018253 
0.014269 
0.009828 
0.005008 
0.000000 


PHASE  ANGLE 


1-53756Q 

1.553782 

1.569994 
-1.555159 
-1.538568 
-1.521979 
-1.505386 
-1.488472 
-1.471694 
-1.454839 
-1.437988 
-1.420986 
-1.404230 
-1.387417 
-1.370529 
-1.353691 
136905 
120334 


-1 

-1.3 
-1.237225 
-1.270993 
-1.116021 
-0.981865 
-0.872103 
-0.786891 
-0.724721 
-0.682756 
-0.660769 
0.273944 


NORM. PRESS 

7.554441 
7.460239 
7.36955" 
28034 


192361 
107576 
024406 
J42234 


7 
7 
7 
7 

6. 

6.863304 
6.785389 
6.709667 
6.635274 
6.56295" 
6.49195 
6.422443 
6.354450 
6.283318 
" ,223366 
,159449 
,096689 
,035670 
,464833 
,917839 
,325723 
.644764 
2.849282 
1.962534 
1.000000 
0.000001 


66 


WEDGE  ANGLE  =    10.00 
SOURCE  DISTANCE=0.}5 
RH01/RH02=  0.80 
K1X   =   19.44 


SOURCE  ANGLE  =    2.50 
RECEIVER  DISTANCE=    2.25 
Cl/C2=  1.10 


SHORE  ,DISTANCE=  0.00 
ALPHA/K2=   0.0001 


REC.POS 

REC. ANGLE 

\ 

0.00 

0.10 
0.20 

3 

4 

0.30 

5 

6 

818 

7 

0.60 

8 

0.70 

9 

0.80 

10 

0.90 

11 

1.00 

12 

1.10 

13 

1.20 

14 

1.30 

15 

1.40 
1.50 

16 

17 

1.60 

18 

1.70 

\l 

1.80 

1.90 

21 

2.00 

22 

3.00 

23 

4.00 
5.00 

If 

6.00 

26 

7.00 

27 

8.00 

28 

9.00 

29 

10.00 

PRES. AMPLITUDE 


Q 

.224816 

0 

.224336 

0 

.223838 

0 

.223305 

0 

•222733 

0 

.222113 
.221472 

0 

0 

.220805 

0 

.220096 
,219349 

0, 

0, 

,218573 
,217751 
,216893 

0. 

0, 

0, 

,216000 

0, 

,215070 

0. 

,214095 

0. 

,213090 

0. 

,212046 

0, 

210963 

0. 

209839 

0. 

208671 
194767 

0. 

0. 

176765 

0, 

154712 

0. 

128893 

0. 

099744 

0. 

068024 

0. 

034481 

0. 

000000 

PHASE  ANGLE 

1.084841 
1.092171 


,099430 

,106651 

113891 

■  121132 

,128275 
,135338 
,142416 
149437 
156374 
163238 
170144 
176949 
183646 
190398 
197003 
203545 
210009 
216443 
222806 
231938 
333861 
77699 


1.413359 
1.441621 
1.461364 
1.473173 
-1.243458 


NORM. PRESS 

6.519929 

6.505989 
6.491551 
6.476104 
459502 
4415! 
42291 
40358! 


6 

6 

6 

6 

6. 

6 

6.338860 

6.315019 

6.290153 

6 

6 

6 

6 

6 

6:08556( 

6.051688 

5.648450 

5.126376 

4.486338 

3.738033 

2.892702 

1.972762 

1.000000 

0.000001 


1.383023 
.361359 


164237 
J37287 
209000 
179859 
149565 
118167 
>5f 
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WEDGE  ANGLE  =    10.1 
SOURCE  DISTANCED. "?! 
RH01/RH02=  0.80 
K1X   =   19.44 


REC.POS 

RE C. ANGLE 

I 

0.00 

2 

0.10 

3 

0.20 

4 

0.30 

s 

0.40 

f 

0.50 
0.60 

8 

0.70 

9 

0.80 

10 

0.90 
1.00 

11 

12 

1.10 

13 

1.20 

if 

1.30 

1.40 

16 

1.50 

17 

1.60 

18 

1.70 

19 

1.80 

20 

1.90 

21 

2.00 

22 

3.00 

23 

4.00 

it 

5.00 

6.00 

26 

7.00 

27 

8.00 

28 

9.00 

29 

10.00 

SOURCE  ANGLE  =    2.50 
RECEIVER  DISTANCE=    1.13 
Cl/C2=  1.10 


SHORE  DISTANCE=  0.00 
ALPHA/K2=   0.0001 


PRES. AMPLITUDE 

1.900229 
1.923264 
1.944952 
1.965253 
1.984233 
2.001945 
2.018349 
2.033438 
2.047371 
2.060028 
2.071478 
2.081736 
2.090869 
2.098818 
2.105642 
2.111385 
2.115994 
2.1195SO 
2.122012 
2.123518 
2.123900 
2.074631 
1.939167 
1.730041 
1.459312 
1.138785 
0.780462 
0.396619 
0.000001 


PHASE  ANGLE 


■0 


0.125102 

0.119954 

114742 

109495 

104208 

-o!088218 

-0.062850 

-0.077489 

-0.072146 

-0.06o788 

-0.06144o 

-0.056122 

-0.050809 

-0.045508 

■0.040234 

~  ,034997 

029773 

024576 

,019431 

,  ,029804 

0.074109 

0.112477 

0.144356 

0.169405 

0.187392 

0.198255 

-1.531829 


NORM. PRESS 

4.791066 
4.849142 
4.903825 
4.955010 
5.002365 
5.047521 
5.088881 
5.126925 
5.162055 
5.193968 
5.222836 
5.248699 
5.271727 

5!308975 
5.323455 
5.335076 
5.344040 
5.350248 
5.354045 
5.355009 
5.230785 
4.889239 
'  ,361967 
679378 


871230 
967786 
000000 
000002 
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WEDGE  ANGLE  =    10.00 
SOURCE  DISTANCED. 75 
RH01/RH02=  0.80 
K1X   =   19.44 


SOURCE  ANGLE  =2.50 
RECEIVER  DISTANCE=    0.56 
Cl/C2=  1.10 


SHORE  DISTANCE^  0.00 
ALPHA/K2=   0.0001 


REC.POS 
1 


REC. ANGLE    PRES .AMPLITUDE    PHASE  ANGLE     NORM. PRESS 


0.00 

5. iq 

0.20 
0.30 
0.40 
0.  50 
0.60 
0.70 
0.80 
0.90 


00 
10 
20 
30 
40 
50 


1.60 


70 

80 

:E8 

,00 

,00 

,00 

00 

00 


8.00 

9.00 

10.00 


3.992624 

4.016742 
4.038875 
4.059057 
4.077303 
4.093628 
4.107439 
4.120035 
4.130817 
4.139806 
4.147007 
4.152471 
.156213 
.158244 
.158599 
.157314 
.154387 
.149846 
.143713 
.136028 
•126792 
s. 953794 
5.647243 
3.224684 
2.703265 
2.100711 
1.435872 
0.728641 
0.000000 


0.306449 
Q. 310540 
0.314667 
0.318822 
0.323001 
0.327198 
0.331284 
0.335505 
0.339736 
0.343969 
0.348202 
0.352429 
0.356651 
0.360864 
0.365063 
0.369249 
0.373416 
0.377565 
0.381693 
0.335796 
0.389876 
0.428945 
0.464082 
0.494513 
0.519774 
0.539597 
0.553331 
0.562396 
-0.847966 


5.479548 

5.512648 

5.543023 

5.570723 

5.595763 

5.618168 

5.637122 

5.654409 

5.669208 

5.681543 

5.691426 

5.698924 

5.704061 

5.706848 

5.707335 

5.705572 

5.701555 

5.695323 

5.686906 

5.676359 

5.663683 

5-426256 

5.005542 

4.425613 

3.710009 

2.883053 

1.970616 

1.000000 

0.000000 
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WEDGE  ANGLE  =    10.00 
SOURCE  DISTANCED.  75 
RH01/RH02=  0.80 
K1X   =   19.44 


tEC.POS 

REC. ANGLE 

\ 

O.OQ 
0.10 

3 

0.20 

8 

0.30 

0.40 
0.50 

8:18 

9 

0.80 

10 

0.90 

11 

1.00 

12 

1.10 

13 

1.20 

it 

1.30 

1.40 

16 

1.50 

17 

1.60 

18 

1.70 

19 

1.80 

20 

1.90 

21 

2.00 
3.00 

22 

23 

4.00 

24 

5.00 

25 

6.00 

26 

7.00 

27 
28 
29 


8.00 
9.00 
10.00 


SOURCE  ANGLE  =    2.50 
RECEIVER  DISTANCE=    0.28 
Cl/C2=  1.10 


PRES. AMPLITUDE 


642356 
636151 


0.629941 
0.623727 
0.617507 
0.611286 
605059 
598825 
592588 
586344 
580099 
573846 
567591 
561331 
555062 
548791 
542518 
536235 
529949 
523657 
517360 
454118 


390374 
155 


326 L  , 

,  261508 

0.196416 

0  " 

0 

0 


13110! 

06560! 
000001 


SHORE  DISTANCE=  0.00 
ALPHA/K2=   0.0001 


PHASE  ANGLE 

-0i270972 
-0.270501 
-0.270035 
-0.269567 
-0.269111 
-0.268658 
-0.268203 
-0.267760 
-0.267313 
"66874 
66441 


-0 
-0 
-O.i   014 


-0.26558! 
-0.265165 
•0.264745 
264339 
263923 
263521 
'63120 


■0 
■0 
■0 
-0 


-0.262721 
"116 


-0.2590L. 
-0.255795 
-0.253081 
-0.250935 
-0.247807 
-0.246445 
-0.245626 
-0.026186 


NORN. PRESS 

9.791314 
9.696724 
9.602075 
9.507351 
9.412549 
9.317716 
9.222797 
9.127778 
9.032702 
8.937539 
8.842341 
8.747032 
8.651687 
8.556263 
8.460711 
8.365122 
8.269497 
8.173734 
8.077916 
7.98200' 
7.88602 
6.922033 
5.950404 
4.971522 
3.986116 
2.99393' 
1.99845. 
000000 


0.000009 


70 


APPENDIX  D 
TABLES 

TABLE  1 

RECEIVER  DISTANCE  AT  THE  FIRST  TRANSITION  POINT, 
FOR  CONSTANT  pj  p2  =  0.80,  Cj  c2=1.10 


P  =  6°.     K,X 

=  32.61 

R,  =  1.10 

R,  =  1.20 

R,  =  1.30 

R,  =  1.40 

R,  =  1.50 

Y  =  P4 

5.2 

5.0 

5.9 

7.3 

9.8 

Y  =  P2 

5.2 

4.7 

5.7 

6.8 

9.9 

Y=3P4 

4.0 

4.5 

5.4 

7.0 

10.5 

P=10°-    K,X=  19.44 

R,=0.80 

R,=0.90 

R,  =  1.00 

R,  =  1.10 

R,  =  1.20 

Y  =  P4 

17.5 

24.0 

33.0 

52.0 

72.0 

y=P/2 

12.5 

17.0 

24.0 

42.0 

60.0 

Y=3P4 

10.6 

17.0 

22.0 

40.0 

58.0 

P=150'   K,X=  12.79 

R|  =  1.30 

R,  =  1.40 

R,  =  1.50 

R,=  1.60 

R,  =  1.70 

Y=p4               S.l 

8.2 

8.7 

9.8 

11.4 

Y  =  P2 

6.9 

7.0 

7.8 

8.9 

10.9 

Y=3P4 

9.8 

12.0 

15.0 

19.0 

29.0 
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TABLE  2 

RECEIVER  DISTANCE  AT  THE  FIRST  TRANSITION  POINT, 
FOR  CONSTANT  Pj  p2  =  0.80,  Cj  c2=  1.20 


P  =  6°,     K,X=22.53 


R,  =  1.30 

R,  =  1.40 

R,=  1.50 

R,  =  1.60 

R,  =  1.70 

Y=p4 

3.37 

3.6 

3.92 

4.4 

4.96 

Y=P2 

3.155 

3.42 

3.78 

4.36 

4.96 

Y=3P4 

3.07 

3.34 

3.74 

4.2 

5.1 

P=100'    K,X=  13.43 


R,=0.80 

R,=0.90 

R,  =  1.00 

R,  =  1.10 

R,  =  1.20 

Y  =  p4 

6.1 

11.3 

18.0 

40.0 

60.0 

Y=p2 

6.5 

8.2 

14.5 

25.0 

40.0 

Y=3p4 

5.05 

5.45 

12.0 

18.0 

30.0 

p=  15°'    K,X  =  8.84 


R,=0.80 

R,=0.90 

R,  =  1.00 

R|  =  1.10 

R,  =  1.20 

Y  =  p4 

36.0 

40.0 

52.0 

60.0 

64.0 

Y  =  P2 

24.0 

26.0 

32.0 

46.0 

54.0 

Y=3p4 

19.50 

23.0 

28.0 

3S.0 

5S.0 
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TABLE  3 

RECEIVER  DISTANCE  AT  THE  FIRST  TRANSITION  POINT, 
FOR  CONSTANT  pj  p2  =  0.90.  Cj  c2=1.10 


P=  6°'   K,X 

=  32.61 

R,  =  1.10 

R,  =  1.20 

R,  =  1.30 

R,  =  1.40 

R,  =  1.50 

Y  =  P4 

3.66 

4.1 

4.9 

6.2 

11.0 

Y  =  P2 

3.35 

3.9 

4.7 

7.0 

no 

Y=3P4 

3.24 

3.76 

4.75 

7.3 

9.8 

P=10°,     K,X=  19.44 

R,  =  l.30 

R,  =  1.40 

R,  =  1.50 

R,  =  1.60 

R,  =  1.70 

Y=p4 

9.8 

11.0 

14.0 

40.0 

60.0 

Y=P2 

7.6 

8.6 

11.0 

19.0 

40.0 

Y=3P4 

6.7 

7.6 

10.0 

21.0 

45.0 

73 


APPENDIX  E 
GRAPHS  OF  Rj  VERSUS  R2  AT  THE  FIRST  TRANSITION  POINT 
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Figure  E.l     Rj  vs  R2  at  the  first  trans.points,  for  p=  6°,  p1/p2-s0.80,  Cj/c2-  1.10 
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Figure  E.2    Rj  vs  R2  at  the  first  trans,  points,  for  p=  10°,  p1/p2  =  0.80,  Cj/^-  1.10 
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Figure  E.3     1^  vs  R2  at  the  first  trans,  points,  for  p=  15°,  p1/p2  =  0.80,  Cj/c2=  1.10 
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Figure  E.4    Rj  vs  R2  at  the  first  trans,  points,  for  P=  6°,  Pj/p2  =  0.80,  Cj/c2  =  1.20 
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Figure  E.5     K{  vs  R2  at  the  first  trans,  points,  for  P«=  10°,  p1/p2  =  0.80,  Cj^-  1.20. 
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Figure  E.6    Rj  vs  R2  at  the  first  trans,  points,  for  p=  15°,  pj/P2  =  0.80,  Cj/C2  =  1.20 
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Figure  E.7     Rj  vs  R2  at  the  first  trans,  points,  for  p=  6°,  p1/p2  =  0.90,  Cj/c2=  1. 
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Figure  E.8     Rj  vs  R2  at  the  first  trans,  points,  for  p=  10°,  p1/p2  =  0.90,  Cj/c2=  1.10 
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